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Abstract - We use the notions of reflexivity and of 
reflexive dimensions in order to introduce probabil- 
ity measures for lattice polytopes and initiate the in- 
vestigation of their statistical properties. Examples 
of applications to discrete geometry include a study 
of randomness of self-duality of reflexive polytopes 
and implications for expectation values of the num- 
bers of such polytopes in higher dimensions. We 
also discuss enumeration problems and related algo- 
rithms and point out interesting open problems. In 
this context we define the notion of IP-confined poly- 
topes. Our new results include the list of IP-simplices 
in 3 dimensions that are not IP-confined. The main 
motivation for the study of these issues comes from 
applications in algebraic geometry and string theory. 

1. Introduction 

Let us consider d-dimensional convex polytopes 
A, which are convex hulls of a finite numbers of points 
in M.'^. Obvious notions of randomness in terms of 
random selections of coordinates of the vertices pre- 
fer simplicial polytopes, whose duals are simple [1]. 
Probability measures based on such concepts hence are 
not invariant under the elementary duality exchanging 
vertices with bounding hyperplanes. 

If we focus on lattice polytopes, whose vertices be- 
long to a fixed lattice Z'' C M'^, we are usually inter- 
ested in affine unimodular equivalence classes, whose 
representatives differ by lattice automorphisms. We 
will need some special classes of such polytopes: 

Definition 1.1. An IP-polytope is a lattice polytope 
A with exactly one interior lattice point, which we can 
choose to be the origin € Z''. 

The polar (or dual) polytope A* of a polytope A C M'' 
with in its interior is defined by 

A*^{yeZ'^ I y-x>-l Vx e A}. (1.1) 

A refiexive polytope is an IP-polytope A for which A* 
is also a lattice polytope. 

Remark 1.2. The number of refiexive polytopes with 
given dimension is finite. Because eq. (1.1) defines an 
involution, the polar polytope A* is refiexive if and 
only if A is refiexive. Refiexivity of an IP polytope is 
equivalent to the property that all facets have lattice 
distance 1 from the origin, i.e. there are no parallel 



lattice hyperplanes between the origin and the sup- 
porting hyperplanes of the facets. 

Definition 1.3. Haase and Melnikov [2] have shown 
that every lattice polytope A is (isomorphic to) a face 
of a refiexive polytope. The reflexive dimensions rd(A) 
is the smallest dimension for which such a refiexive 
polytope exists. 

Since the number of refiexive polytopes is finite 
we can define probability measures on the set of all 
polytopes that are uniform in fixed dimension and 
thus respect the polar duality of refiexive polytopes 
exchanging vertices and facets. Similarly, we can use 
the fact that the number of polytopes A with fixed 
refiexive dimension rd(A) is finite to define measures 
on the set of all lattice polytopes. The numbers of 
polytopes with rd( A) < 4 are given in table 2 below. 

The content of this paper is organized as follows. 
In section 2 we summarize results concerning the enu- 
meration of refiexive polytopes and discuss some rel- 
evant algorithms. In section 3 we extend the polar- 
ity of refiexive polytopes to an involutive duality of 
IP-confined polytopes in cZ > 4 dimensions and dis- 
cuss various enumeration problems. Section 4 contains 
statistical considerations and observations, which are 
used in an attempt at estimating the number of refiex- 
ive polytopes in 5 dimensions. 

Acknowledgements. The author would like to ac- 
knowledge discussions with Benjamin Nill and support 
by the Austrian Research Fund FWF under grant no. 
P18679-N16. 

2. Reflexivity and weight vectors 

The classification of refiexive polytopes in two di- 
mensions can be done by hand on a piece of paper 
with the result shown in fig. 1. The interest in ex- 
plicit enumerations in higher dimensions is mainly mo- 
tivated by the relation between lattice polytopes and 
algebraic geometry, and in particular by the result of 
Batyrev [3] that generic hypersurfaces in toric vari- 
eties are Calabi-Yau if and only if the toric ambient 
space is defined in terms of a refiexive polytope. On 
the algebro-geometric side of this correspondence the 
polar duality corresponds to an exchange of complex 
structure and Kahler moduli with important applica- 
tions to enumerative geometry [4,5]. In terms of string 
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Figure 1. The 16 reflexive polygons in 2 dimensions: 3 are maximal and the last 4 are selfdual. 



theory, this amounts to determining instanton correc- 
tions from a duality called mirror symmetry, which ex- 
changes matter with anti-matter in particle physics [5] . 

The algorithm that has been used for the com- 
plete enumeration of reflexive polytopes in 3 and 4 
dimensions [6, 7] is based on the idea to determine 
the maximal objects in terms of the polar minimal 
ones [8] and on an algorithm for the enumeration of 
all possible linear relations among vertices of the min- 
imal polytopes [9]. For the case of two dimensions 
the three maximal reflexive polytopes are drawn, to- 
gether with the dual minimal ones, in the top row of 
fig. 1 and it is easy to check that all others are sub- 
polytopes thereof. If the maximal polytope is sim- 
plicial then the coefficients in the linear relation be- 
tween the vertices of the dual minimal lattice poly- 
tope, which are chosen positive integral without com- 
mon divisor, are called weight system or weight vec- 
tor w — {wq, . . . ,Wd)- For the two simplicial mini- 
mal polytopes in fig. 1 the weight vectors are (1,1,1) 
and (1,1,2) and they can be obtained as Wi = D qi 
where D is the common denominator of the barycen- 
tric coordinates Qi of the origin with respect to the 
vertices of the minimal simplices, i.e. ^ QjVj — and 
~ 1- The polar maximal polytopes can be re- 
covered as Newton polytopes (i.e. exponent vectors of 
the monomials) of generic quasi-homogeneous polyno- 
mials of degree D w.r.t. the weights Wi. The 10 points 
of the first maximal polytope, for example, correspond 
to the 10 monomials of degree 3 in 3 variables. 

Non-simplicial minimal IP-polytopes are convex 
hulls of a collection of lowcr-dimcnsional minimal IP- 
simplices and thus define combined weight systems, or 
weight matrices, whose rows consist of weight vectors 
augmented by zeros for the vertices that do not be- 
long to the respective simplex. For the third minimal 
polytope in fig. 1, i.e. the square, the matrix of linear 
relations reads (q ^ g ^) , while a minimal pyramid in 
three dimensions is described by a weight matrix of 
the form C"" ""^^ ""l,' ° ° ). 

\W21 U'24 



Definition 2.4. For wj £ N with gcd{wj) = 1 let 
A(w) = ConvHull({m}) denote the Newton polytope 
of a generic quasi-homogeneous polynomial J2m. '-™ ^J^' 
of degree D = J2j regarded as a lattice polytope 
on the intersection of the lattice U^'^^ with the hyper- 
plane m-w = D. Since Wj > and rrij > the Newton 
polytope A(w) has at most one interior point, namely 
rrij = 1. We call w an IP-weight-vector if A(w) is an 
IP-polytope. 

Remark 2.5. It follows already from a more general 
result of Lagarias and Zicglcr [10] that the number 
of IP polytopes, and hence the number of IP weight 
vectors, is finite in fixed dimension. A constructive 
proof has been found in terms of an efficient algorithm 
for the enumeration of all IP- weight- vectors [9, 11]. 

Remark 2.6. Since all IP-polygons are reflexive, lin- 
ear relations of vertices of IP-simplices in two dimen- 
sions are always IP- weight- vectors. For d > 3, on the 
other hand, the latter form a proper subset of the set 
of all linear relations among simplicial IP-polytopes. 
The lowest degree examples of non-IP-weights of IP 
simplices in 3, 4 and 5 dimensions are Pi.5^6.8[20], 
Pi,i,5,5,8[20] and P2,2,2,3,3,ii[23], respectively. Table 1 
contains the complete list for c? 3 in boldface. Note 
that Wj < D/2 for all linear relations of IP-simplicies 
because for wq > wi + . . . + Wd the simplex (uj) with 
'^WjVj — has at least two interior lattice points, 
namely and —vq — '^^j^j with cq — ^jf- — 1 and 
Cj = -jf-Vj for j > 0, so that — 1 and Cj > 0; if 
Wq = D/2 then —vq is in the interior of a facet. 



Theorem 2.7 (Skarke [9]). The Newton polytope /S.{w) 
of an IP weight vector w — {wq^ . . . ,Wd) is always 
reflexive in d < A dimensions. More generally, for 
weight-matrices Wij with l<i</, l<j<J and 
Di = X^ Wij the Newton polytopes /S.{W) of generic 
polynomials of multi- degrees Di are always reflexive 
whenenver W is an IP-weight-matrix (i.e. when A{W) 
is an IP polytope) if the dimension d~ J — I of A.{W) 
obeys d < 4. 

The first counterexample \s w — (1,1,1,1,1,2) with 
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Table 1. The 104 IP-simplices in d = 3 correspond to 95 transversal and 9 (boldface) non- IP- weights. 



degree D = 7, whose 5-dimensional Newton polytope 

A{w) has one facet with lattice-distance 2 from the 
interior point and hence is not reflexive. 

Remark 2.8. The number of IP weight vectors is 
Nip{l) = 1, Nip{2) = 3, iV/p(3) = 95 and Nip{4) = 
184026 for d < 4. By extrapolation we expect that an 
enumeration might be possible in d = 5 but certainly 
not above because the numbers become too large. For 
enumeration problems discussed in the next section it 
would hence be important to find an efficient general- 
ization of Skarke's algorithm [9] for IP weight vectors 
with a bounded number of lattice points in their defin- 
ing minimal simplex. 

Definition 2.9. A weight vector w is called transver- 
sal if a generic quasi-homogeneous polynomial /a(«,) = 
SmGA Cm Y[j °^ degree D is transversal df{x) = 
=> X = 0, i.e. the gradient only vanishes at the 
origin xj = 0. Due to the Euler formula f{x) = 
^QjXj df /dxj with qj = Wj/D the equation / = 
thus defines an isolated singularity [12]. 

Remeirk 2.10. Transversal weight vectors are a sub- 
class of the IP weight vectors which concide with the 
reflexive weights for d < 3 and imply reflexivity for d = 
4 [9]. They have been enumerated for d < 5 [13-16] 
with Nt{1) = 1, ^t(2) = 3, 7Vt(3) = 95, 7Vt(4) = 
7555 and A/t(5) = 1100 055. The reflexive weights 
are much more numerous than the transversal ones 
already in 4 dimension. Nevertheless, only 252 933 of 
the transversal weights in 5 dimensions, about 23%, 
are reflexive. The above example (1, 1, 1, 1, 1, 2) of a 
non-reflexive IP weight is also the simplest example of 
a non-reflexive transversal weight. 

The complete algorithm for the enumeration of re- 
flexive polytopcs is compiled in [11]. A somewhat sub- 
tle detail is the possible choice of a sublattice on the 
side of the maximal polytope, which is constrained by 
the fact that the surviving lattice points need to span 
an IP polytope on the sublattice. An implementation 
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Table 2. Numbers of lattice poly- 
topes with reflexive dimension < 4. 



of the algorithm and of additional applications to com- 
binatorics and algebraic geometry has been published 
in the program package PALP [17]. With some fur- 
ther extensions, which are implemeted in the current 
version 1.1, this package has been used for all com- 
putations of results that are report in this note. In 
particular, we enumerated all lattice polytopes with 
reflexive dimension^ rd(A) < 4, whose numbers are 
listed in table 2. 

3. IP-confined polytopes 

The algorithm for the enumeration of reflexive 

polytopes is simplified in d < 4 by the fact that New- 
ton polytopes A(H^) of IP weight matrices are auto- 
matically reflexive so that all maximal reflexive poly- 
topes are of this form. Since this is no longer true for 
d > 5 we are lead to consider a slightly larger class of 
IP polytopes. 

Definition 3.11. Let A := ConvHull(A*nZ'^) be the 

convex hull of the lattice points in A* . An IP polytope 
is called IP-confined (IPC) if A is a lattice polytope. 

^ The program and many data are available at [18]. Basic de- 
scriptions of extended /experimental features are available with 
the help option -x. In particular, facets arc enumerated with 
poly.x -F and class. x -A creates lists of affinc normal froms. 
The third column of tabic 2 can be obtained with 

$ echo "41111 /Z2: 1 1" > 3d. in; cws.x -w3 » 3d. in 
$ cws.x -c3 >> 3d. in; class. x -po 3d.ref 3d. in 
$ class. X -b -pi 3d.ref I poly.x -fF I class, x -f -A -po 3rd. 2d 
$ class. X -B -pi 3rd. 2d I poly.x -fF I class, x -f -A -po 3rd. Id 
after fetching and compiling PALP, e.g. as in (4.6) below. 



The polytope A c A is called the IPC-closure of A 
and A is called IPC-closed if A = A. 

Remark 3.12. The operation A — > A extends the 
polar duality of reflexive polytopes to an involution 
on the set of IPC-closed polytopes. (Reflexive poly- 
topes are obviously IPC-closed; according to Skarke's 
theorem [9] the reverse is true for d < 4.) 

Corollary 3.13. Each IP- confined polytope (and hence 
also each reflexive polytope) is a subpolytope of a max- 
imal IPC-closed polytope Am- These maximal poly- 
topes are of the form Am = ConvHull(A(VF) n F) 
where W is an IP-weight-matrix with rank(W^) < d 
and V (Z is a sublattice of the natural lattice of 
the Newton polytope for which A{W) fl T still has an 
interior point. 

Proof. The proof is analogous to the one for the reflex- 
ive case [11] and uses the generalized duality A i— > A 
of definition 3.11. The polytopes A that are dual to 
the maximal IPC-closed polytopes are minimal in the 
set of IPC-closed polytopes, but not necessarily mini- 
mal as IP polytopes. For every minimal IP polytope 
A' C Am, however. Am is the IPC-closure, because 
otherwise Am would not be maximal. A weight ma- 
trix W with Am = ConvIIull(A(VF) n F) can hence 
be constructed by regarding A' as the convex hull 
of lower- dimensional IP-simplices Sj (which are IP- 
confined since the duals to the restrictions to the linear 
subspaces of Sj are projections of Am, which provide 
confining IP polytopes for the simplices Ej). Each 
simplex in the decomposition can be chosen such that 
it contains at least two vertices of A' that are not 
contained in a different simplex of the decomposition. 
The maximal rank d hence occurs if the maximal poly- 
tope is a hypercube. □ 

The numbers of reflexive polytopes in d dimen- 
sions are the diagonal entries of table 2. A double- 
exponential ansatz N^i « 2^ ~^ due to Skarke pro- 
vides a good fit to these data and would predict « 
1.2 • 10^^ and TVe ~ 2.1 • 10^^ in 5 and 6 dimensions, 
respectively. If this gives the correct order of magni- 
tude an enumeration would already be hopeless in 5 
dimensions because the result could not be stored on 
any existing medium (for the 4d-case we needed, on 
average, about 20 byte, or five 32-bit integers, to store 
the data of a reflexive pair). Bounds on the numbers 
of lattice points have been obtained in [19] . 

While the case of 4-dimcnsional reflexive poly- 
topes is of direct relevance to algebraic geometry and 
string theory because the corresponding toric hyper- 
surfacc^s are Calabi-Yau 3-folds [3] , higher-dimensional 
cases are important because F-theory applications [20] 
in physics require (elliptically fibered) Calabi-Yau 4- 
folds, and, more generally, reflexive polytopes in n + r 
dimensions are related to complete intersection n-folds 



of codimension r [21 23]. From a pragmatical point 
of view, in both contexts the polytopes with a small 
number of points are the ones that we are most inter- 
ested in because they encode the combinatorial data 
for manifolds with (moderately) small Picard num- 
bers. This suggests a different setting for enumera- 
tion attempts where we ask for the lists of reflexive 
polytopes with a certain (maximal) number of points 
rather than a fixed dimension. Instead of descending 
from maximal objects, which is not feasible for d > 4, 
one might hence add points to a list of minimal ob- 
jects. 

Question 3.14. Is there an efficient algorithm for 

the enumeration of weight vectors w that define IP- 
confined simplices J2j '^j^j = with a bounded num- 
ber of lattice points? 

While we expect the number of IP-weight- vectors 
to be too large for d > 5 for a complete enumeration, 
we also expect that the enumeration of all reflexive 
polytopes with less than, say, 12 or 15 lattice points 
should be feasible if the questions 3.14 has a positive 
answer. Here the natural objects of interest are IP- 
confined lattice polytopes because every IPC-closed 
polytope, and hence every reflexive polytope, can be 
constructed by successively adding lattice points to 
the minimal objects in a finite number of ways as de- 
termined by the relevant weight vectors. 

4. Statistics and numerical experiments 

We now turn to statistical considerations that may 
provide additional information in cases where a com- 
plete enumeration is not possible. The first question 
we want to address is whether the property of self- 
duality of reflexive lattice polytopes is random in the 
sense that the number of selfdual objects can be esti- 
mated by a probabilistic calculation. 

Let N be the number of reflexive polytopes in 
flxed dimension d. Then the number of (formal) dual- 
ity assignments with S self-mirrors, i.e. the number of 
involutions with S flxed points, is (^) for the choice 
of S self-mirrors times {N - S -1)-{N - S -^)- . . .-i-l 
for the possible selections of the remaining dual pairs, 
hence 



ns 



S\-2 



^(^)! 



with S-N ^21. (4.2) 



Let Zf^ = '^g^^p^ns be the total munber of involu- 
tions. The asymptotic expansion for large N can be 
derived from the generating function [24, section 3.8] 



- 2^ N\ 

N>0 



(4.3) 



Since it happens Z^-i times for the Zn involutive 

permutations of N objects that a given object is a 
fixed point, we obtain the following formula for the 
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Figure 2. Predictions for reflexive transversal Newton polytopes with < 36 . . . < 65 points. 



expectation value [25] 



(4.4) 



Assuming a uniform probability distribution of invo- 
lutions we thus expect 



1 1 

3 sjN 



N I 



(4.5) 



selfdual polytopes for large iV. As shown in table 3 
this roughly explains the size but underestimates the 
correct numbers for c? < 4. 



d 


1 


2 


3 


4 


Nnid) 


1 


16 


4319 


473800776 


A^'sclfdual 


1 


4 


79 


41710 




1 


3.6 


65.2 


21766.5 



Table 3. Numbers of selfdual poly- 
topes and probabilistic expectations. 

As this result seems to indicate that a statisti- 
cal approach makes sense we now want to use similar 
considerations for predicting the number of reflexive 
polytopes in 5 dimensions on the basis of incomplete 
lists. In a random set of p > ^/N reflexive polytopes 
in fixed dimension the formula (5*) ~ \fN implies that 
we expect s « ^/N -p/N self- mirror polytopes. This 
leads to the prediction N w [p/ s)"^ if we find s self- 
dual polytopes in the sample. Similarly, if we ignore 
the relatively small number of self-mirrors for large N 
and count the number m of mirrors pairs in a sample 
of p polytopes we obtain the prediction N « p^ /{2m). 
If we increase the size of a random sample we expect s 
to grow linearly and m quadratically with p, and more 
precisely s w \/2m. 

Actually, the formula N w p^ / {2m) has been used 
already serval years ago when we enumerated the re- 
flexive polytopes in 4 dimensions [7], which required 
two years of program improvements and computation 
time after the 3-dimensional case [6]. We thus could 
check the sufficiency of the implemented data struc- 
tures for the storage of the result at an early stage of 
the project. The starting point of the calculation was 
the list of 308 weight matrices (206 weight vectors and 



102 matrices with 2 < rank < 4) of maximal reflex- 
ive polytopes Aa/ for which the reflexive subpolytopes 
were computed in the order of an increasing number 
of lattice points. The first Newton polytope in this list 
is defined by the weight vector (3, 3, 4, 4, 10) with de- 
gree 24. It has 47 lattice points and 6 vertices. After 
fetching and compiling PALP its subpolytopes can be 
computed with the following commands. 



# directory for binary files (check $PATH) 

# working directory 



$ Bdir-$HOHE/bin 
$ Wdir=$PWD 
$ cd /tmp 

$ wget hep . itp . tuwien, ac . at/~kreuzer/CY/palp/palp-i . 1 
$ gunzip palp-* , tar .gz : tar -xvf palp-*, tar; cd palp; 
$ make ; mv * , x $Bdir ; cd $Wdir 
$ 

$ echo "24 3 3 4 4 10" I class, x -f -po /tmp/zbin,47 



tar,gz 



# fetch 

# unpack 
# and compile PALP 



800kR-1658 IIHB 4215klP 1042kNF-23k 9_46 Vl7rl7 f28r27 85b24 120s 120u 25n 
24 3 3 4 4 10 R-798878 +1658sl hit=0 1P=4215623 NF-1042005 (23261) 
,,, 1181m+14s 9851469b u7 pp/2m-2,68615e+08 pp/ss=3 , 26616e+09 

(4.6) 

Within two minutes we thus obtain the values 269 mil- 
lion for p'^l2m with p = 798878 and to=1181 and 3.25 
billion for p^ / with s = 14 on a standard 3GHz PC 
(cf. the last output line) and hence a good approxima- 
tion of the correct value 473 800 776 with a production 
rate of about 7000 polytopes per second (back in 1998 
the CPU time was almost 1 hour). 

For 5 dimensions it is, of course, much harder 
to get a reliable statistics. With an expectation of 
N « 10^^ according to Skarke's guess, the storage of 
\/iV polytopes would already required some 30GB of 
disk space so that we can only go above that value 
by 1-2 orders of magnitude with currently available 
hardware. Nevertheless, it should be possible to ei- 
ther verify that the number is not much smaller or to 
get a reasonable prediction if it is. For a first attempt 
we defined data samples in terms of the Newton poly- 
topes of transversal reflexive weight vectors^ ordered 
according to increasing numbers of lattice points. For 
the increasing series data samples consisting of all re- 
flexive subpolytopes of transversal Newton polytopes 
with < 36 . . . < 65 points the predictions are plotted 
in fig. 2 

The result is obviously inconclusive, but certainly 
compatible with the guess 10^^. Unfortunately the 



cf, http: //quark, itp. tuwien. ac . at/~ kreuzer/W/5dTR* . 



bias of the a priori independent predictions /2m and 
/s^ is strongly correlated in our data samples. For 
the 9.025 • 10^ polytopes of the largest sample with 
m = 86323 and s = 354, which occupies 239 GB of 
disk space, we find p^ /2m « 4.7 • 10^^ and p^/s^ « 
6.5-10". 

5. Conclusions 

In this note we determined all lattice polytopes 
with reflexive dimension rd < 4 and discussed enu- 
meration problems and algorithmic aspects with ap- 
plications to algebraic geometry and string theory. We 
pointed out the need for an efficient algorithm for the 
enumeration of IP weight vectors w with a bounded 
number of lattice points in the convex hull of the sim- 
plex defined by the linear relations ^ WjVj = 0. Such 
an algorithm could be used for the enumeration of re- 
fiexive polytopes with fixed number of points rather 
than fixed dimension. 

We introduced the concept of IP-confined poly- 
topes, which are a subclass of IP polytopes, and ex- 
tended the polar duality of reflexive polytopes to IPG- 
closed polytopes. Maximal IPG-closed polytopes con- 
tain all reflexive polytopes in arbitrary dimensions and 
hence lead to a simplification of the classification pro- 
gram. In turn, we pointed out the existence of IP- 
simplices that are not IP-confined and enumerated 
them for the case of 3 dimensions. A constructive 
classification of such simplices in higher dimensions is 
another interesting open problem. 

We suggested a statistical approach to the enu- 
meration of reflexive polytopes which should at least 
allow us to obtain probabilistic lower bounds, depend- 
ing essentially on the size of the available hard-disks 
for storage of the data. As a first attempt we con- 
structed a data-base containing about 9-10^ pairs of 
refiexive 5-dimensional polytopes, which can also be 
used to produce incomplete lists of polytopes with re- 
flexive dimension 5. Studies of correlations of polytope 
data like /-vectors (numbers of faces) can thus be ini- 
tiated and may be useful for selecting appropriate data 
samples for statistical applications. 
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